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ABSTRACT 

On UT 2009 January 16, we observed a white light megaflare on the dM4.5e 
star YZ CMi as part of a long-term spectroscopic flare-monitoring campaign to 
constrain the spectral shape of optical flare continuum emission. Simultaneous 
[/-band photometric and 3350-9260A spectroscopic observations were obtained 
during 1.3 hours of the flare decay. The event persisted for more than 7 hours 
and at flare peak, the [/-band flux was almost 6 magnitudes brighter than in the 
quiescent state. The properties of this flare mark it as one of the most energetic 
and longest-lasting white light flares ever to be observed on an isolated low-mass 
star. We present the [/-band flare energetics and a flare continuum analysis. 
For the first time, we show convincingly with spectra that the shape of the blue 
continuum from 3350A to 4800A can be represented as a sum of two components: 
a Balmer continuum as predicted by the Allred et al radiative hydrodynamic flare 
models and a T ~10,000K blackbody emission component as suggested by many 
previous studies of the broadband colors and spectral distributions of flares. The 
areal coverage of the Balmer continuum and blackbody emission regions vary 
during the flare decay, with the Balmer continuum emitting region always being 
significantly (~3-16 times) larger. These data will provide critical constraints for 
understanding the physics underlying the mysterious blue continuum radiation 
in stellar flares. 

Subject headings: stars: flare — stars: late-type — stars: atmospheres — stars: 
individual (YZ CMi) 
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Introduction 



Due to their strong and persistent surface magnetic fields, active M dwarfs can produce 



frequent flares, sometimes lasting for many hours and re aching 
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uminosities that approach 



20091 ) . During both solar and stellar flares, emission is seen as line and contin 



uum radiation from X-ray to radio wavelengths, with a prominent mode of radiative energy 
release occurring in the blue and near-ultraviolet (NUV) continuum, which is commonly 
referred to as the white light continuum. Unfortunately, the spectral components which 
comprise the white light continuum are not well-constrained, and therefore the physics of 
the underlying flare mechanism that gives rise to such a conspicuous signature remains a 
mystery. 

Significant effort has gone into characterizing the white light conti nuum in stella r flares , 
beginning with a single component (hydrogen recombination) model (IKunkell Il969l Il970l ). 
A two-component spectral model consistin g of hyd r ogen recombination and an impulsively 
heated photosphere was first proposed by IKunkell (119701 ) as a possible e xplanation of the 
sprea d in broadband flare colors. Time-resolv ed multi-channel photometry (IMochnacki fc Zirin 
19801 ) and spectrometry ( jKahler et al.lll982l ) were obtained longward of 4200 A for flares on 
YZ CMi and UV Ceti; at peak, these flare spectra were fit with a Planck function of T = 
7400— 9500K, thereby showing evidence of a hot optically thick component present during the 
impulsive phases of flares. Instrumental effects shortward of the B aimer jump (3646A) pre- 
vented a precise characterization of the Balmer continuum in the data of IMochnacki fc Zirin 
Jl980h . 



Higher reso lution spectra in the blue /N UV were obtained during larg e flares on AD Leo 
(3560 A- 4440A; lliawlev fc Pettersenl [l~99~l and CN Leo (3050A- 3860A; iFuhrmeister et al. 
20081 ) showing a nearly continuous rise into the NUV with aT~ 8400 — 11, 300K blackbody; 
surprisingly no discontinuity at the Balmer jum p was seen in either of these events (see also 
Eason et all 1 19921 ; iGarcfa- Alvarez et all 120021 ). The white light continuum has therefore 
been attributed to a single dominant blackbody component with T ~ 9000 — 10, 000K (e.g., 
Hawley et al. l boo3h . In contrast, flare models predict spectra that are dominated by hydrogen 



bound-free continu a, which provide a poor match to the observations (IHawley &; Fisherlll992 
Alfred et al.lbood . hereafter A06). 



A multi-component model was reintroduced by IZhilyaev et al.l (120071 ). who used high- 
cadence UBVRI photometric observations of the dM3.5e star EV Lac to suggest that the 
white light emission consists primarily of blackbody radiation at flare peak and hydrogen 
continuum during the flare decay. However, significant ambiguity results when using broad- 
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band photometry to characterize the white light continuum. A06 showed that a flare model 
spectrum that included a large Balmer jump, when convolved with broadband UV + UBVR 
filters, actually resembles the shape of a T ~9,000K blackbody. To unambiguously identify 
the spectral components present during stellar flares, we have therefore begun to compile 
a large catalog of time- resolved optical/NUV flare spectra on several nearby flare stars. In 
this Letter, we present time- resolved spectra obtained during a large flare on the dM4.5e 
star YZ CMi and use results from the recent radiative hydrodynamic (RHD) flare models of 
A06 to show that both Balmer continuum and hot blackbody components are present in the 
white light continuum emission. An extensive paper describing our flare catalog, including 
discussion of the flare emission lines and comparison with models will be forthcoming. 



2. Observations 

2.1. NMSU 1-m Photometry 

Photometric observations were obtained in the Johnson U band on UT 2009 January 
16 with the New Mexico State University 1-m Telescope at the Apache Point Observatory 
(APO). The 1-m is operated robotically a nd is open for observa tions when the ARC 3.5-m 



Telescope, located at the same site, is open ( IHoltzman et al.ll2010l ). The dM4.5e star YZ CMi 



was observed for nearly 8 hours with an exposure time of 10s; with readout the photometric 
cadence ranged between 20 - 27s. We performed differential photometry on YZ CMi using 
the comparison star HD 62525. The night was clear and photometric based on the data from 
nearby comparison stars. The [/-band light curve of YZ CMi is shown in Figure []Ja. 



2.2. APO 3.5-m Spectroscopy 

We obtained NUV/optical spectroscopy for 1.3 hours during the decay phase of the flare, 
denoted by vertical lines in Figure []Ja, using the Dual Imaging Spectrograph (DIS) with the 
ARC 3.5-m Telescope at APO. We used the B400+R300 gratings, providing wavelength 
coverage from 3350A— 9260A. The observations were taken through a 1.5 " slit which yielded 
resolutions of R ~ 625 (at 4000A) and R ~ 980 (at 6563A), determined from the FWHM 
of the arc lamp lines. The slit was oriented at the parallactic angle in order to compensate 



for atmospheric differential refraction (see iFifippenkol 119821 ) . We observed YZ CMi with 



exposure times of 10s, and with readout the total cadence was ~28s. In total, 163 flare 
spectra were obtained simultaneously in both the blue and red bandpasses. The data were 
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reduced using standard IRAF0 procedures, wavelength calibrated using HeNeArHg lamp 
spectra, and flux-calibrated by compa ring observ ations of the spectrophotometric standard 
star Feige 34 to the accepted fluxes of lOkel (ll990l ). To correct for variable grey slit loss (e.g., 
due to guiding errors), we also scaled each spectrum during the flare using the [/-band light 
curve. A spectrum of YZ CMi in its quiescent (non-flaring) state was obtained with the same 
instrumental setting on UT 2008 November 24 and was used for comparison (see §3.2) J^| 



3. Analysis 

3.1. [/-Band Flare Energetics 

This flare event is remarkable in its total energy, peak luminosity, total duration, an d 
emission morphology. A small precursor, sometimes seen prior to large flares ( IMoffettlll974j ). 
began at 04:14:54UTC which increased the [/-band flux by a factor of ~2, and represents the 
onset of the flare event. At 04:32:00UTC, there was a ~10 min impulsive rise to flare peak, 
where the [/-band flux was a factor of ~200 times the quiescent value. The decay phase 
included many smaller events, and the [/-band flux was still elevated by over a magnitude 
above quiescence at the end of the night. In §4, we discuss the extremely long duration of 
this event. 

The total [/-band flare e nergy is used to quantify the energy release in continuum 
radiation during stellar flares ( Moffctt|l l974f ). We integrate under the light curve to obtain 



an equivalent duration ( IGershberg 



19721 ) of 4.25 x 10 5 s, which is the amount of time that the 



star would spend emitting at the pre-flare level to produce the same total amount of energy 
as the flare. Multipl ying the equiv alent duration by the non-flaring [/-band luminosity of 
4.00 xlO 28 ergs s" 1 flMoffettl Il974f ) gives E UJlare > 1.7 x 10 34 ergs (this is a lower limit 
since flare emission is still seen at the time the light curve ends). At peak, the flare was 
emitting with a [/-band luminosity of 8.3 xlO 3 ergs s _1 , or ~37% the stellar bolometric 



luminosity (Mb i= 10.25, iReid &: Hawleyll2005l ). This flare was ~100 times more energetic 
and persisted for more than 10 times as long as the largest [/-band flare eve nt observed on 
this star during 55 hours of monitoring from the large statistical flare study of lMoffettl (Il974f ) 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation 

2 The quiesc e nt sp ectrum was convolved with the Johnson B and V bands (using fil ter curves from 
Maiz Apellaniz ( 20061) ) and scaled to match the flux from the literature (V — 11.2, B = 12.8; Reid fc Hawlev 
2005) 
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and lLacy et al.l (119761) . A l thoug h it is unknown if the power law fit to the flare-frequency 
distribution of lLacy et al.l ( 119761 ) holds at high flare energies, an extrapolation predicts that 
flares with [/-band energies > 1.7 x 10 34 ergs occur at a rate of ~once per month. 



3.2. The Flare Continuum 

The spectral data were obtained while the [/-band emission was still elevated 15-37 
times (3-4 mag) over quiescence. Complex morphology is seen in the [/-band light curve 
during this time period (see Figure [Ih), allowing us to study the flare continuum variations 
during several smaller impulsive events. Figure [Tb shows a spectrum (t = 76.6min after the 
flare start) from 3350A— 5500A near the beginning of our spectroscopic observations during 
the [/-band decay from a large secondary emission peak. A quiescent spectrum of YZ CMi 
(see §2.2) is also shown in Figure [Tb and has been subtracted from the flare spectrum to 
give a "flare-only" spectrum. Our spectra have the highest time resolution with high signal 
to noise in the blue/NUV (S/N > 50 at 3600A) ever obtained during a large stellar flare. 
They are also unique in having simultaneous high quality photometric information so that 
we may relate the spectral evolution to the light curve morphology. 

We have for the first time convincingly identified two distinct continuum components 
that simultaneously contribute to the optical/NUV flux in flare spectra. The first component 
dominates the spectrum from 4000-4800A and exhibits the rising trend of a T ~10,000K 
blackbody, which is plotted on the spectrum in Fig ure Hb using a fillin g factor, X#b, of 



0.22% of the projected stellar disk area. We follow lHawley et al.l ( 120031 ) to constrain the 



filling factor of blackbody emission as a function of time using the equation: 

FflM = X BB (t)^irB x {Tfl) ( 
where Fn \ is the flare flux observed at Earth in the continuum windows fr om 3995-4020A, 



4140-4210A, 4425-4450A, and 4600-4800A; R* = 0.3 x R Sun , d = 5.97 pc (IReid fc Hawlev 



2005|) and T fl = 10, 000K. 



From Figure [Tb, it is clear that the blackbody fit does not explain all of the continuum 
emission at wavelengths shorter than ~3800A. Starting near H10 (3799A), there is a rise to 
H15+H16, the last two apparent Balmer lines in emission. Blueward of ~3700A, the flare 
spectrum becomes flat into the NUV, and there is excess flux visible above the 10,000K 
blackbody at these wavelengths. 



A06 presented RHD flare simulations for a medium (F10) and a strong (Fll) impulsive 
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flux of nonthermal electrons into an M dwarf atmosphere. The A06 models predict strong 
Balmer continuum radiation with a noticeable Balmer jump at 3646A, together with a pho- 
tosphere that is warmed by ~1200K to a temperature of 4600K. We isolated the Balmer 
continuum component of the A06 Fll model spectrum taken at the last time-step of their 
simulation (t = 15.9s), and fit it to the observed flare spectra. Figure [lb shows the flare 
spectrum from Figure [lb with the 10,000K blackbody component subtracted and the model 
Balmer continuum component scaled to the excess emission in the 3600A < A < 3646A 
region. The shape of the observed excess emission is well- described by the Fll flare model 
prediction (shown in red). Although there is no abrupt discontinuity indicative of a Balmer 
jump in our observed spe ctrum, we posit that blending (possibly due to Stark broadening; 



Donati-Falchi et al.lll985l ) of the higher order Balmer lines in the 3646A < A < 3800A region 



effectively hides the sharp discontinuity that is seen in the model spectrum. 



3.3. Relative Filling Factors 

The filling factor for the Fll model spectrum is calculated by scaling the model flux in 
the 3600-3646A region {Fx^odei — 1-14 x 10 7 ergs cm 2 s 1 A 1 ) to the flux excess observed 
at Earth after subtracting the 10,000K blackbody. The contribution function for the Balmer 
continuum in the Fll model shows that the emission arises from a relatively thin layer in 
the atmosphere (J. Alfred, private communication); it is therefore reasonable to approximate 
changes in the filling factor to changes in the emitting area of the Balmer continuum. 

In the upper panel of Figure [Tbl (black points), we show the ratio of the filling factors 
for a 10,000K blackbody and the model Balmer continuum component. The error bars are 
obtained by using 9000 and 11,000K temperatures for the blackbody. The ratio for the 
sample spectrum in Figure [lb is X BB /X Mode i = ClOig^. The filling factor ratio indicates 
that the two emission components exhibit significantly different time evolution, with the 
blackbody source emitting from a region ~3-16 times smaller than the source of the Balmer 
continuum. 

The inferred area filling factor (Xbb) of the blackbody component is shown separately 
(red points) and exhibits very similar morphology to the [/-band light curve (reproduced 
in the lower panel). Although we expect the total continuum flux to follow the [/-band 
since we scaled the spectra to achieve absolute flux calibration (§3.2), it is surprising that 
only the blackbody component follows the U bandit It appears that the deviations from 
the general flare decay seen in the U band primarily arise from the changing area of the 



3 These trends are also visible in the data before the flux calibration scaling was applied. 
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blackbody emission region. When the deviations are strongest near the peaks at t ~ 95 
min and t ~ 130 min, the spectra are increasingly dominated by the b l ackbody component. 



This i s cons istent with the flar e -peak spectra from lKahler et al.l (11982). lHawlev fc Pettersen 



( 1199 it ), and iFuhrmeister et al.l (120081 ) and the photometry of IZhilyaev et al.l (120071 ). We also 
note preliminary evidence of an apparent increase in the best-fit blackbody temperature to 
T ~ 13, 000K at the peak near t ~ 130 min. 

The derived filling factor ratios should be regarded as illustrative since our Balmer 
continuum spectrum from the Fll model of A06 assumes that the nonthermal electron flux 
is constant at 10 11 ergs cm -2 s _1 . Also, we use the model flare spectrum from only one 
time-step in the RHD simulation. It is likely that the nonthermal electron flux varied both 
with time and over the geometry of the flaring region. Any further assessment is beyond the 
scope of this Letter, but we note that the less energetic A06 F10 flare spectrum at t — 230.0s 
gives nearly 10 times larger filling factors for the Balmer continuum. 



3.4. Balmer Line and Continuum Radiation 

The time-evolution of the Balmer continuum closely follows the hydrogen Balmer lines. 
We show the variation in the Balmer continuum (using the integrated flux in a 30A window 
centered at 3615A) compared to the H7 line flux and the U band in the lower panel of 
Figure [TH. The fluxes are normalized to the last spectral observation we obtained. The 
Balmer continuum sh ows a slow decay through out our observations, as is typical of the 



Balmer line radiation (lHawlev fc Pettersenlll99lt ). This overall decay is probably extended 



emission from the previous, much more luminous, impulsive events seen in the [/-band light 
curve. The larger variations in the Balmer continuum compared to the H7 emission may be 
partially due to small errors in isolating the Balmer continuum by subtracting the underlying 
blackbody continuum (the H7 flux is well measured as it has a local continuum). As noted 
above in the filling factor discussion and evident by comparing the upper and lower panels 
of Figure [TH, the Balmer line and continuum radiation exhibit quite different time evolution 
compared to the blackbody emission component. We propose two possible explanations: 

(a) The hydrogen-emitting regions exhibit a delayed and more gradual response to each suc- 
cessive impulsive (blackbody) event. In this scenario, the gradual increase in hydrogen 
Balmer emission around t = 122 and t = 137 min are associated with the impulsive 
events in blackbody emission (upper panel) at t = 117 min and t = 127 min, respec- 
tively. Also, Balmer emission around t = 80 min and t = 105 min is elevated yet only 
slowly decreasing, and could be gradually declining emission from the large [/-band 
peaks at t ~ 60 min and t ~ 95 min, respectively. 
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(b) The time evolution of the blackbody emission and the Balmer emission may be anti- 
correlated. If true, this may provide an important constraint on the origin of the 
blackbody emission. The possible anti-correlation is most apparent during the ?7-band 
rise starting at t — 124 min. In ~ 5 min the Balmer continuum decreases by ~40% 
while the blackbody filling factor almost doubles. Hydrogen continuu m dimming is 
predicted by flare models (lAbbett fc Hawleyl Il999t lAllred et al.l 120051 ). but only for 
very short durations (0.1s) at the onset of a flare. 



More detailed models will be required to understand these complex variations, and in 
particular, the seemingly anti-correlated trends between the blackbody and Balmer contin- 
uum and lines during the time intervals when the blackbody filling factor is increasing in 
Figure HH. 



4. Summary and Speculation 



We observed a white light megaflare on the dM4.5e star YZ CMi in the U band and with 
simultaneous optical/NUV spectroscopy. The [/-band energetics and light curve morphology 
qualify this flare as an extraordinary a n d rare event, simil ar to the ~ 6.6 x 10 34 erg flare on 



YZ CMi seen bv lAndrews et al.l (119691 ). lLovelll (1196a ). and|Kunke| ( 119691 ). Following a solar 
analogy, we speculate that the flaring region on YZ CMi was a complex arcade of sequentially 
reconnecting magnetic loops. Each reconnecting loop accelerated a beam of nonthermal 
electrons that impacted the lower atmosphere, producing the observed blue/NUV line and 
continuum emission. The sum of a large number of individual emitting regions may have 
enabled this flare to persist for such an unusually long time. 

Using high time resolution, high signal-to- noise spectra, we have shown that the blue/NUV 
flare continuum radiation can be explained as a sum of a T ~10,000K blackbody compo- 
nent and a Balmer continuum component, with only the blackbody emission present from 
4000-4800A. The relative filling factors of the two components indicate that the Balmer 
continuum comes from a larger region, plausibly originating from the flaring loops at chro- 
mospheric heights, where the Balmer lines originate. This is consistent with the height of 
the Balmer continuum emission derived from solar flare data ([Hudson et al.ll201Cf ). and also 
with the height of formation predicted by the RHD model of A06. 

We found that the blackbody emission arises from a region ~3-16 times smaller in 
area than the Balmer continuum emission region. This blackbody emission may possibly 
originate in concentrated magnetic footpoint regions in t he lower atmosphere , similar to the 



localized areas that emit in white light during solar flares (IMetcalf et al.ll2003l ; iFletcher et al. 
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20071 ; llsobe et al.ll2007t I Jess et al.ll2008l ). which are often spatially and temporally coincident 
with impulsive heating by nonthermal electrons inferred from hard X-ray observatio ns (e.g., 



Rust fe Hegwerlll975l : lHudson et al]ll992 iNeidig fc Kandll993l : Fletcher et al.ll2007l ). In ac- 



cordance with this scenario, a large complex of photospheric hot spots may have been created 
during the first impulsive events in the YZ CMi ?7-band light curve. As the nonthermal elec- 
tron beams weakened (in energy and/or in spatial extent) so did the areal coverage of these 
hot spots. This gradual decay extends into the time covered by our spectral observations 
and is seen in the overall decline of the blackbody areal coverage in the upper panel of Figure 
[TH. During our spectral observations, a few new hot spots may have been created at the 
footpoints of newly reconnected loops, causing the transient increases seen in the [/-band 
and in the effective area of the blackbody emitting region. 

We see strong evidence through all of our observations for a blackbody continuum 
emission component with an approximate temperature of 10, 000K. The persistence of the 
hot blackbody emission indicates that a continual source of particle acceleration and plasma 
heating likely still exists during the decay phase of the flare. Yet, the A06 models, which 
employ a nonthermal electron beam as is seen on the Sun, predict the photosphere of an 
M dwarf to be heated by at most ~1200K during a large flare. The physical mechanism 
which generates the strong blackbody emission (presumably from a more strongly heated 
photosphere) therefore remains unknown. The possible anti-correlation between the Balmer 
continuum and blackbody emission components may provide an important clue to the nature 
of the heating mechanism that is responsible for the blackbody emission. 
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Fig. 1. — (a) The flare [/-band light curve with vertical lines indicating the period of spec- 
troscopic observations, (b) A flare spectrum at t — 76.6min after flare start. The quiescent 
spectrum is shown in purple for comparison. The continuum from 4000 - 4800A was fit with 
a T = 10, 000K blackbody and a filling factor of 0.22%. There is excess emission above 
the blackbody blueward of 3800A. (c) The 10,000K blackbody component has been sub- 
tracted from the spectrum in panel (b) and the Fll flare spectrum from A06 fits the excess 
(Balmer) continuum at wavelengths shorter than 3646A. (d) (upper panel) The filling factor 
ratio (X BB /X Mode i, black points) and the inferred area coverage of the blackbody component 
(red points), (lower panel) The time-evolution of the [/-band (left axis), and the Balmer 
continuum and H7 fluxes (right axis). 



